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SUMMARY  Setting  up  and  solving  ballistic  type  rocket  kinetic 
equations  in  inertial  coordinate  systems  with  directions  which  do  not 

vary,  it  is  possible  not  to  consider  dragging  inertial  or  centrifugal 
forces  or  coriolis  inertial  forces  and  set  up  inertial  guidance  on  a 
convenient  and  simple,  quick  and  clear  foundation.  This  is  different 
from  what  has  been  published  in  the  past,  where  studies  have  been  done 
of  various  types  of  flight  mechanics  problems  in  coordinate  systems 
which  turn,  following  along  with  the  earth.  We  illucidate  the 
initial  take  off  conditions  as  rockets  leave  the  surface  of  the  earth, 
and  analyze,  in  detail,  the  air  or  atmospheric  effects  associatea  with 
air  being  pulled  along  by  the  gravity  of  the  earth  as  it  turns  and  the 
treatment  methods  for  aerodynamic  forces  on  rockets  moving  in  an 
inertial  space.  l*ic  derive  a  formula  for  calculating  the  geographical 
location  of  movement  parameters,  used  for  rockets  in  inertial 
coordinate  systems,  and  a  formula  for  dynamic  coordinate  system 
movement  parameters. 

KEY  TERMS  Inertial  Guidance,  Coordinates,  Flight  Mechanics, 
Kinetic  Equation  .  - 


As  far  as  the  literature^^'^'^'^^and  publications  associated 
with  the  kinetics  and  kinematics  of  domestic  and  foreign  studies  of 
lifting  rockets  and  other  similar  types  of  spacecraft  are  concerned, 
although  all  start  out  through  classical  mechanics,  in  the  end, 
however,  they  all  infer  the  handling  of  aerodynamic  forces  in 
coordinate  systems  which  turn  together  along  with  the  earth  as  well  as 
various  similar  types  of  forces  and  motion  parameters.  Moreover,  it 
is  necessary  to  consider  centrifugal  inertial  forces  and  coriolis 
inertial  forces  as  well  as  other  similar  types  of  virtual  forces. 
Speaking  in  terms  of  the  researchers  who  have  handled  guidance  and 
navigation  systems  for  rockets  and  other  similar  types  of  spacecraft, 
if  one  analyzes  motion  characteristics  of  inertial  instruments,  it 
becomes  even  more  complicated.  In  many  years  of  research  and 
practical  application,  we  took  rocket  kinetics  problems  and  solved 
them  in  an  unmoving  inertial  space,  obtaining  satisfactory  results. 


I.  SEVERAL  COORDINATE  SYSTEMS 
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1.  The  Earth  Centered  Coordinate  System  E,X  , Y  ,Z  . 

^  ^ 

The  origin  point  E  is  placed  at  the  geometrical  center  of  the 

earth.  The  EX  axis  in  the  plane  of  the  equator,  from  the  center  of 
o 

the  earth  E,  points  toward  the  point  of  intersection  of  the  prime 
meridian  and  the  equator  and  is  positive  in  direction.  The  EY^ 
axis,  in  the  plane  of  the  equator,  is  perpendicular  to  the  EX^  axis 
and  is  positive  in  direction  toward  the  outside.  The  EZ^  axis, 
along  the  axis  of  the  earth's  autorotation,  forms  a  clockwise 
orthogonal  coordinate  system  with  EX  and  EY  .  This  is  as  shown 
in  Fig.l. 

2.  Guidance  Coordinate  System  0^,X^,Y^,Z^ 

The  origin  point  O  is  the  launch  point  of  the  rocket  on  the 

earth's  surface.  The  0  X  axis,  within  the  horizontal  plane  that 

o  o  , 

passes  through  0^,  points  toward  the  direction  of  the  flight  and  is 

positive.  The  included  angle  or  angle  of  offset  between  it  and  the ' 

due  north  direction  of  the  meridian  group  tangent  is  the  azimuth 

angle.  The'O^Y^  axis  follows  along  a  plumbline  upward  and  is 

positive.  The  axis,  in  the  horizontal  plane  that  passes 

through  0^,  forms  a  clockwise  orthogonal  coordinate  system  with 

0  X  and  0  Y  .  At  the  instant  that  the  rocket  leaves  the 
o  o  o  o 

launch  pad  (time  t=0sec),  one  takes  the  direction  O^X^Y^Z^  is 
pointing  in  inertial  space  and  firmly  fixes  it  invariable.  It  then 
forms  an  inertial  coordinate  system.  See  Fig.l. 


3.  Launch  Coordinate  System  or  Dynamic  Coordinate  System  OXYZ. 

At  the  instant  the  rocket  takes  off  (t=0®),  the  origin  points  0 

and  0  ,  the  OX  axis  and  OX,  the  OY  axis  and  O  Y  ,  and  the 
o  o  o  o  o 

OZ  axis  and  OqZ^,  are  all  superimposed  on  each  other.  After  the 
rocket  leaves  the  ground  and  takes  off,  OXYZ  is  firmly  connected  to 
the  surface  of  the  earth.  Moreover,  it  turns  together  with  the  earth. 
See  Fig.l. 


4.  Spacecraft  Coordinate  System 

The  origin  point  Oj^  is  positioned  at  the  center  of  mass  of  the 
rocket  body.  The  O^Xj^  axis  runs  along  the  main  axis  of  inertia  of 
the  rocket,  points  toward  the  nose,  and  is  posi_tive.  0^2 in  the 
secondary  plane,  is  perpendicular  to  the  plane  of 

symmetry.  This  is  as  shown  in  Fig. 2. 


v. 


Fig.l  The  Relationships  of  Such  Coordinate  Systems  as  EX^Y^Z^, 
and  0XY2 

O'  O  O  -• 


Fig. 2  The  Relationship  of  and 
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5.  Airflow  Coordinate  System  02^2^2^! 

The  origin  point  O2  is  conincident  with  the  center  of  mass  of 
the  rocket  body  Oj^.  The  axis  02^2  line  with  the  air  flow 

velocity  vector  corresponding  to  the  rockets  The  02^2 
in  the  main  rocket  plane  goes  upward  and  is  positive. 

The  02^2  forms  a  clockwise  rectangular  coordinate  system  with 

the  02^2  *^2^2  takes  the  O2X2Y2Z2 

coordinate  system  and  rotates  it  about  the  O2Y2  axis  through  the 
lateral  slide  angle  ,  and,  again,  rotates  it  about  the 
axis  through  the  angle  of  attack  a,  then  it  is  coincident  with 
0,X,Y,Z.  See  Fxg«3. 


Fig. 3  The  Relationships  of  0^,X.,Y.,Z.  and 
02,X2,,Y2,22  .  1  1  i  i 


The  unit  vectors  on  the  axses  of  each  coordinate  system  are  taken 

as  a  unit  1,  adding  the  corresponding  coordinate  system's  lower  corner 

subscripts  to  express  that.  For  example,  1^„,1^  ,  and  1  are 

,  ox  oy  oz 

the  unit  vectors  of  the  guidance  or  inertial  coordinate  system 

Y_Z_  on  the  three  axses  0  X  ,  O  Y  ,  and  0  Z  . 

0000  000  o  o  o 


II.  KINETIC  EQUATIONS 


Taking  out  the  influences  of  other  celestial  bodies,  one  only 
considers  the  effects  of  earth's  gravity,  G  ,  The  high  speed  flow 
of  combustion  gases  sprayed  out  by  the  rocket  engines  produces  as  an 
effect  the  thrust  vector  on  the  rocket.  Its  magnitude  is  related 
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to  the  amount  of  mass  consumed  as  fuel  in  a  unit  time  G=dG/dt  (kg/s), 
the  engine  jet  tube  characteristics,  and  the  height  h  at  which  the 
rocket  flies  off  the  surface.  Normally,  this  can  be  expressed 


as 


[1,2,3] 


(2.1) 


On  the  basis  of  the  principles  of  momentum,  the  rocket's  instantaneous 
acceleration  vector  in  the  inertial  coordinate  system 
OqX^Y^Z^,  V^®dV ^/dt,  satisfies  the  equation  below 


(2.2) 


Force  vector  F  is  the  vector  sum  of  thrust  F_,,  control  force  F  , 

2  ^ 

and  the  pneumatic  force  F  .  g  =9.80665  (m/s  )  and  is  the  mass 

<3  e 

conversion  constant.  is  the  engine's  specific  impulse. 

^e(unclear)  engine  combustion  gas  flow  aperture  area.  P/Pq 

is  the  ratio  of  the  gas  pressufe  and  the  surface  atm.ospheric  pressure 

P^,  which  varies  with  the  height  h  from  which  rockets  fly  off  the 

surface.  If"  one  takes  a  certain  vector,  such  as  F,  and  uses  a  line 

matrix  to  express  it,  its  symbols  and  the  vector  symbols  are  the  same. 

T 

One  then  has  F=(FjjFyF2]  .  The  upper  right  superscript  T 
expresses  transposed  positions.  F^^,  F^ ,  and  F^  are  the 
components  of  F  on  the  and  axses.  When 

using  the  acceleration  of  gravity  g=  comparative 

force  or  perveived  acceleration  W=(Wj^WyW2  3  to  express 
equation  (2.2),  one  obtains 


(2.3) 


In  this 


Ip 

.^F. 


The  force  vector 


F“Fr+F*+Ff 


(2.3.1) 


32 


The  thrust  force  and  control  force  are  normally  given  in 
cubic  coordinate  systems.  Moreover/  the  aerodynamic  force  Q2  in  the 
gas  flow  coordinate  system  02X2^2^2  can  be  transformed  into 
^l^l^l^l  expressed.  In  this  way/  the  comparative  force 
components  on  the  main  inertial  Oj^X^^  /  O^^l'  °1^1 

axses  are  respectively 

^Xl™'^iFtTX'^FiKX+Qlx) 

^  (y.d.li 


Wri’^’^iFtrr+FtKr+Otr) 


(2.4.2) 


Wxt^^-^iFtrx+FtKi+Qtz) 


The  comparative  force  vector  The  .  '  . 

relationship  between  the  unit  vectors  1q  and  Ij^  of 
O^X^Y^Zo  and.  can  be  used  to  link  up  or 

associate  the.  orthogonal  transformation  matrix  D/  which  is  nothing 
else  than 


l.-D-l, 


(2.5) 


Because  of  thiS/  the  comparative  force  vector  W  in  the  guidance 
coordinate  system  O^X^Y^Z^  is  capable  of  being  aided  by 
using  the  comparative  force  vector  Wj^  in  to  express 

as 


W^D-W, 


(2.6) 


to 


Gauging  or  measurement  assemblies  used  in  the  guidance  and  stability 
systems  of  such  types  of  carrier  or  delivery  craft  as  rockets  are,  it 
goes  without  saying,  gyroscope  stabilized  platforms  or  sensitively 
connected  inertial  assemblies.  The  various  elements  of  the  orthogonal 
transformation  matrix  D,  dj^^  (i  =  l,2,3,  j=l,2,3),  are  all  capable  of 
being  arrived  at.  For  example,  in  Fig. 2,  rotating  the  pitch  angle  cp 
around  the  axis  0q2q,  and  again,  rotating  the  yaw  angle  ^around 
the  transitional  axis  OY'j^,  and,  finally,  revolving  the  roll  angle 'V 
around  the  axis  in  that  case,  d^^^  is  then  possible  to  make 

use  of  to  obtain  the  sine  and  cosine  functions  of  Euler  angles  <p  , ' 
and  V  . 

Taking  rockets,  when  they  act  as  rigid  bodies,  the  Euler 

equations  for  the  revolved  center  of  mass  0.  are  normally  set  up  in 

[61  ^ 
a  cubic  coordinate  system^ 


JxAtx  +  Uz-Jr)caiz&tr’*Mix  (2.7.1) 

Jr^xr-¥{J  x'-J  i)o%xCitt'^M%r  (2.7.2) 

/jjwix +  (/r  ~/r)oir<»ix“Mix  (2.7.3) 

In  the  equations,  Jj^,  Jy,  and  are  the  instantaneous  momenta 
of  rotation  of  the  rocket  about  the  main  axses  of  momentum 
and  Mj^y,  ^re  component  moments  of 

force  along  the  three  axses  discussed  above.  c*j 

are  components  of  angular  acceleration  for  rotation  of  the  rockec 
about  Oj^Yj^,  and  The  angular  acceleration 

vector  is 


Oi-*Coijr  Oir  (2.8) 

This  article  only  illucidates  movements  of  the  center  of  mass.  It 
does  not  discuss  the  kinetic  equations  for  motions  around  center  of 
mass . 
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III.  INITIAL  CONDITIONS  AND  AERODYMANIC  FORCES 


In  the  inertial  coordinate  system  OXYZ^,  as  far  as 

o  o  o  o 

study  of  the  flight  kinetics  of  rockets  and  other  similar  carrier  or 
delivery  vehicles  is  concerned,  the  most  important  thing  is  to 
illucidate  the  methods  of  handling  the  effects  of  aerodynamic  forces 
with  which  the  atmosphere  acts  on  rockets  in  high  speed  flight  as  it 
rotates  along  with  the  earth,  and,  on  the  foundation  of  discussions 
with  fellow  workers,  in  conjunction  with  the  carrying  out  of 
theoretical  analysis,  to  obtain  correct  solutions. 

Observed  in  inertial  space,  rockets,  before  they  take  off,  rotate 
together  along  with  the  earth.  They  are  dragged  along  by  the  force  of 

t 

the  earth's  gravity  and  connected  firmly  to  the  surface  of  the  earth. 
However,  they  already  possess  a  drawing  or  centrifugal  linear 
acceleration.  This  velocity  vector  is  the  initial  velocity 

taking  off  from  earth  under  the  effect  of  the  thrust  of  the  engines. 
The  magnitude  of  is  determined  by  the  vertical  distance  of  the 

launch  point  from  the  earth's  axis  of  autorotation.  In  the  same 
way,  the  atmosphere,  which  wraps  around  the  earth,  is  also  pulled 
firmly  by  the  force  of  the  earth's  gravity  and  follows  the  earth  in 
its  rotation.  Observed  from  inertial  space,  this  atmosphere  also 
possesses  similar  linear  acceleration. 

Looking  from  the  revolving  earth,  the  rocket,  before  it  takes 
off,  is  not  moving.  The  atmosphere  is  stationary.  The  rocket's 
movement  velocity  relative  to  the  atmosphere  is  zero.  As  a  result  of 
this,  there  are  no  aerodymanic  or  pneumatic  force  effects  on  the 
stationary  rocket.  '  This  is  saying  nothing  else  than  that,  if,  when 
one  calculates  the  magnitude  of  the  effects  of  aerodynamic  or 
pneumatic  forces  on  the  rocket  in  inertial  space,  it  is  necessary  to 
take  the  dragging  or  centrifugal  velocity  associated  with  the  rocket 
rotating  along  with  the  earth  and  the  velocity  of  the  atmosphere 
moving  together  with  the  earth  and  subtract  them  from  each  other,  one 
obtains  the  rocket's  velocity  relative  to  atmospheric  movements  as 
zero,  making  the  aerodymanic  forces  borne  by  the  rocket,  when  it  is 
erect  on  the  launch  pad  and  has  not  taken  off  yet,  to  be  zero.  If  one 
takes  the  initial  absolute  velocity  of  the  rocket  before  it 
takes  off,  as  observed  in  inertial  space,  and  takes  it  as  the  movement 
velocity  relative  to  the  atmoshere,  calculating  out  the  effects  of 


8 


pneumatic  forces  on  the  rocket,  in  that  case,  the  rocket  erected  on 
the  launch  pad  will  just  be  pushed  over  onto  the  ground.  However, 
this  does  not  correspond  to  the  facts. 

After  the  rocket,  under  the  effects  of  thrust  F^,  leaves  the 
launch  pad  and  takes  off,  the  absolute  velocity  of  the  motion  in 
inertial  space  gets  greater  and  greater.  The  height  of  the  separation 
from  the  earth's  surface  is  also  able  to  continuously  increase.  The 
magnitude  of  the  aerodynamic  forces  that  it  exerts  at  a  given  point  P 
separated  from  the  center  of  the  earth  E  is  also  determined  by  the 
velocity  of  the  atmosphere  corresponding  to  its  being  pulled  around  by 
the  earth.  It  is  only  in  the  inertial  space  that,  3S 

far  as  the  rocket's  velocity  relative  to  the  atmosphere  is  concerned, 
it  is  appropriate  to  use  this  relative  air  flow  velocity  to  calculate 
the  effects  of  pneumenatic  or  aerodynmaic  forces  on  the  rocket. 

The  components  of  the  earth's  atuorotational  velocity  vector  uJ> 

on  the  various  individual  axses  of  the  inertial  or  guidance 

coordinate  system  O^X^Y^Z^  are: 

■'  o  o  o  o 

.  .  i3x  =  o).'Cos5.-cos‘^,  =  cBl,6x  '  (3.1.1) 


Or  =  (i}t'  sLn5o  =  c}.6r 


(3.1.2) 


Qg  =  (j},-(-cosB,-  sinA,)=<a,-bz 


(3.1.3) 


In  these,  is  the  geographical  latitude  of  the  rocket  launch  point 
Oq.  Aq  is  the  azimuth  angle  of  Oq^q*  O^X^Y^Z^, 

the  absolute  velocity  of  the  air  at  the  place  P(R)  where  the  rocket  is 
located  is  nothing  else  than  the  dragging  or  centrifugal  velocity 
V^.  It  can  be  expressed  as 


V,  =  <a,  X  K  = 


Kz* 

Qx  Oy  S2z 

Rx  Ry  Rz 


(3.2) 


? 


0 


The  magnitude  of  the  radius  vector  from  the  point  P  toward  the  center 
of  the  earth,  R=[R^RyR2]^,  is 

/?=V  (A'.+/?.x)‘  +  (K.+^.r)’  +  (Z.+/?.,)‘  (3.3) 


T 

In  this,  ^OY '  ^OZ  ^  radius  vector  toward  the 

center  of  the  earth  for  the  initial  location  of  the  rocket  Oq  in 
OqXqYqZq.  From  the  origin  point  Oq,  the  position  vector 
pointing  toward  P  is  ^  =[XqYqZqJ^. 

Given  initial  conditions,  that  is  to  say,  the  instant  of  rocket 
take  off  is  t=0sec  and  the  location  relative  to  the  OqXqYqZq 
system  is  zero,  that  is,  Xq=Yq=Zq=0  and  R=R»,  the  initial 
velocity  is  nothing  else  than  the  dragging  or  centrifugal 

velocity  of  point  Oq,  that  is, 

(3.4) 

In  this  way,  the  variable  coefficient  differential  equation  (2.3)  is 
then  capable  of  using  numerical  value  methods  for  its  solution.  The 
absolute  vetocity  of  trie  rocket  Vq  at  a  certain  instant  as  well  as 
the  position  ^  and  R  can  all  be  obtained. 

The  rocket's  velocity  relative  to  the  motion  of  the  earth  is  also 
nothing  else  than  the  motion  velocity  relative  to  the  atmosphere 
and  is  the  difference  between  the  absolute  velocity  Vq  and  the 
dragging  or  centrifugal  velocity  V^(unclear) 


(3.5) 


The  aerodynamic  drag  force  X^  and  the  relative  velocity  are 
mutually  opposite  in  direction,  that  is,  they  are  mutually  opposite  to 
the  direction  of  02^2*  areodynamic  or  pneumatic  lift  (normal 

pneumatic  lift)  Y  is  positive  in  the  direction  along  The 

lateral  pneumatic  force  Z^  is  opposite  in  direction  to  02^2' 

They  are  respectively  used  in  the  several  forms  set  out  below  for 
calculations 

(0 


—Cx'^'Sc 


(3.6) 


Kf  »«CJ  'g'Sc'a 


(3.7) 


Z,~-C*'q-Sc'0 


(3.8) 


In  this,  C^f  C“,  and  are ,  respectively,  the 

aerodymanic  or  pneumatic  drag  coefficient  of  the  rocket,  the 
aerodynamic  lift  coefficient,  and  the  aerodynamic  or  pneumatic 
laterally  directed  force  coefficient.  As  far  as  impact  pressure  is 
concerned,  one  uses  the  form 


(3.9) 


to  calculate  it.  The  air  density  p  changes  along  with  changes  in  the 
height  h  of- the  rocket  as  it  flies  off  the  surface.  is  the 

reference  cross  section  surface  area. 

From  Fig, 3  one  can  see  that  the  airflow  coordinate  system 
9°®®  through  a  rotation  through  the  lateral  slide 
angle  ^  and  the  angle  of  attack  a  to  arrive  at  The 

transformation  between  the  unit  vectors  of  the  two  coordinate  systems 
uses  the  matrix  E  and  is  expressed  as 


(3.10) 


Obviously,  E  is  a  third  order  orthogonal  matrix.  Its  various  elements 
are  functions  of  a  and  /3  : 


~€tt 

pCOs/9*casa 

sina 

—  sin/3*cosa-| 

E- 

- 

—  cos/9*  sina 

cosa 

sin/9‘  sina 

^•1 

^>1 

^sin/9 

0 

008/9  -* 

U 


The  angle  of  attack  a  and  the  lateral  slide  angle  ^  are  calculated  by 
the  use  of  the  formulae  set  out  below; 


(3.12) 


(3.13) 


In  the  cubic  coordinate  system  the  relative  velocity 

V,  =(V,yV,  „V,  is  obtained  by  V,  in  the  O^X^y^Z 
i  Ia  ii  iz  d  o  o  o  o 

system  going  through  transformation.  It  is  capable  of  being  expressed 


as 


V.-D'^V^ 


(3.13) 


Obviously/  are  invariant.  As  far  as  the  aerodynamic  or 

pneumatic  force  airflow  coordinate 

system  02^2^2 ^2  concerned t  it  is  possible  to  make  use  of 

the  matrix  E  <to  transform  it  into  the  cubic  coordinate  system 

^l^l^l^l'  i®' 


O.-E-O. 


(3.14) 


Control  force  vectors  are  always  given  in  the 
^l^l^l^l  coordinate  system.  Practically  speaking,  their 
handling  is  different  because  of  different  types  of  servomechanisms. 
Up  to  now,  we  have  already,  in  O^X^Y^Z^,  completely  obtained 
the  thrust  force  the  aerodynamic  or  pneumatic  force  ,  and 

the  control  force  Fj^.  It  is  possible,  from  equations  (2.4)  and 
(2.6),  to  calculate  out  the  comparative  forces  for  changes  over  time. 


W  and  W 


t(unclear ) * 


When  one  takes  the  earth  to  be  acting  as  a  revolving,  symmetrical 

ellipsoid  sphere,  it  is  possible  to  only  give  consideration  to  the 

gravitaional  force  potential's  second  degree  harmonic  coefficient 

[  5]  ^ 

quantity  .  The  components  of  the  gravitational  acceleration  g, 

which  the  rocket  receives  pointing  toward  the  vicinity  of  the  center 
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of  the  earth  and  which  are  in  directions  along  0  X  ,  0  Y  ,  and 

o  o  o  o 

0^2^  are; 
o  o 


j.- flc •  4’* 0  0 

(y^x.va) 


(3.15) 


In  thiS/  the  constant  of  gravitational  force  is 


Qc'Rc^fM « 


(3.16) 


Rq  is  the  earth's  equatorial  radius,  g^  is  the  acceleration  of 
gravity  on  the  equator.  The  vertical  distance  that  the  rocket's  real 
time  position  P(R)  is  off  the  plane  of  the  equator  is 


C“6x(.^r»  +  ^«») +6r*(K« +/?«r) +6*  *.(Z» +/?*/)  (3.20). 

Th^  point  P's  geocentric  latitude  <p^  is  calculated  from  the 
equation 


Vx- 


(3.21) 


2b^sinCP^  in  g^  does  not  exist  in  OXYZ  when  one  is  studying  the 
motion  equations  for  the  center  of  mass  of  rockets. 

The  application  of  the  various  expressions  set  out  above  is  a 
condition  for  the  complete  solution  of  the  kinetic  equation  (2.3).  It 
is  then  possible  to  precisely  determine,  in  changing  situations,  such 
motion  parameters  as  the  absolute  velocity  and  the  position  R  or  ^ 

IV.  KINETIC  EQUATIONS  FOR  THE  UNPOWERED  SECTION 

As  far  as  the  rocket  or  its  released  nose  cone  is  concerned,  in 
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the  section  of  the  flight  after  the  propulsion  system  is  shut  off  and 
there  are  no  thrust  force  effects  or  attitude  angle  controls,  it  is 
possible  to  take  them  to  be  acting  as  a  point  mass  or  material  point 
particle.  It  is  not  necessary  to  consider  the  angle  of  attack  a  or 
the  lateral  slide  angle  *  The  only  aerodynamic  force  there  is  is 
drag.  Moreover,  it  is  opposite  in  direction  to  , 


Si 


(4.1) 


The  comparative  force  is 


Jt. 

2G, 


-c.rp-yi 


(4.2) 


In  this,  Sj.,  and  are,  respectively,  the  nose  cone's 

aerodynamic  drag  coefficient,  reference  cross  section  surface  area, 
and  its  mass.  Taking  Vl^  and  analyzing  it  along  the  three  axses  of 
the  inertial  coordinate  system  0  X  Y  ,  one' then  has  kinetic 
equations  similar  to  (2.3); 


y _ Ur 

r,y - vv^  '^3y 


(4.3.1) 


(4.3.2) 


(4.3.3) 


In  this,  the  form  expressing  the  acceleration  of  gravity 
T 

g,.g-l  the  same  as  (3.15).  The  initial  conditions  for 

equation  (4.3)  are  the  motion  parameters  V^(tf)  and  §  (tf)  for  the 
instant  when  the  thrust  stops  or  the  instant  tf  when  the  nose  cone  and 
the  rocket  separate. 


In  situations  in  which  the  flight  of  the  nose  cone  reentering  the 

atmosphere  has  attitude  angle  controls,  the  aerodynamic  force  (^2 

contains  the  drag  force  X  ,  the  lift  force  Y  ,  and  the  lateral 

q  q 

force  Zg.  It  is  necessary  to  make  use  of  the  handling  methods  of 
3.  Moreover,  one  must  do  transformations  into  the  coordinate  system 
(for  example,  O^X^Y^Z^)  firmly  connected  to  the  "horse's  head" 
cone.  After  solving  Euler  equations  similar  to  (2.7),  one  obtains  the 
transformation  matrix  D.  In  this  way,  the  aerodynamic  or  pneumatic 
relative  forces  are  also  capable  of  being  accurately  handled  in 


0  X  Y  Z^. 
0000 


V.  MOTION  PARAMETERS  OF  THE  REVOLVING  EARTH 


The  preceding  several  sections  studied  methods  of  handling 
aerodynamic  forces  in  inertial  space  and  the  problems  of  solving 
kinetic  equations.  However,  are  the  results  in  line  with  solutions  of 
kinetic  equations  in  the  dynamic  coordinate  system  OXYZ?.  In  using 

T  t- 

transform  methods,  from  V  =[V  V  V  ]  and  t  =  ' 

,p  o  ox  oy  oz  5 

[X^Y^Z^]  ,  on'e  obtains  dynamic  coordinate  system  motion 
parameters,  and,  in  conjunction  with  that,  after  one  makes  a 
comparison,  one,  then,  has  naturally  dispelled  misgivings.  Refering  to 

Fig.l,  going  through  the  transformations  described  below,  it  is 

possible  to  search  out  the  relationship  between  OXYZ  and 

0000 

OXYZ. 

The  first  step  is  to  rotate  counterclockwise  in 

a  positive  direction  around  O^Y^  an  azimuth  angle  A^  forming 
^0^0  *^0^0'*  tangent  line  O^N  associated  with  the 

meridian  line  of  '  and  the  point  O^  is  duplicated  and  points 

north.  within  the  horizontal  plane  at  point  0  ,  is 

perpendicular  to  the  meridian  plane  and  points  east. 

The  second  step  is  to  rotate  ^qX^'Y^Z^'  counterclockwise 
around  O^Z^'  an  angle  equal  to  the  latitude  B  forming 
^0^0  '  '  ^o  '  ^o  '  *  ^0^0'  parallel  to  the  plane  of  the 

equator,  ^q^o '  '  parallel  to  the  earth's  axis  of  autorotation 
EZ  . 


The  third  step  is,  in  a  situation  in  which  the  three  axses  of 

OX  ''Y  'Z  '  are  maintained  in  an  unchanging  direction,  to 
O  O  O  O  - 

take  and  translate  it  to  the  center  of  the  earth  L  forming 

EX^  "  Y^  '  . 

O  O  O 

The  fourth  step  is  to  take  EX  ''Y^'Z^'  and  rotate  it 
^  O  o  o 

counterclockwise  through  an  angle  around  EX^ ' '  ,  causing 

X  ''EY  '  to  rotate  into  the  meridian  plane  of  point  0. 
o  o  • 


The  fifth,  sixth,  and  seventh  steps  are,  respectively,  the 

inverse  processes  of  the  first,  second,  and  third  steps.  Finally, 

then  takes  O  X  Y  Z  and  transforms  it  to  a  position 
o  o  o  o 

duplicating  OXYZ.  Moreover,  one  obtains  the  transformation  matrix 

It  has  already  been  demonstrated  that  the  matrix  B  is  orthogonal. 

-1  .  .  .  T 

inverse  matrix  B  and  the  transformed  position  matrix  B  are 

equal , 


t),. 

Oi» 

B  = 

b„ 

btt 

&» 

bu 

btf . 

(5.1) 


This  is  capable  of  being  written  to  become  the  equation  below: 


B«B,4-coew,f(/-  B.)  +  sin<y.»-B, 


(5.2) 


In  the  equation,  I  is 
also  both  third  order 


a  third  order 
matrices: 


unit  matrix. 


'6.6. 

bjby 

bJ), 

byby 

.bj>t 

byb. 

'  0 

-b. 

b. 

0 

l_~by 

b. 

B^  and  Bg 

bj}, 
b,b. 
btb.  . 

br 

-b. 

0 


are 


(5.3) 


lb 


one 

B, 

The 


(5.4) 


The  relative  velocity  V  in  the  dynamic  coordinate  system  OXY2  is 
calculated  from  the  equation 


V-BVrf 


(5.5) 


In  this,  Yd* ^'clx'^dy'^dz ^  relative  velocity  expressed 

by  the^use  of  (3.5)  in  O^X^Y^Z^.  The  relative  location 

■IXY2J  in  OXYZ  is  also  capable  of  being  expressed  by  the  use  of 
matrix  B  as 


The  magnitude  r  of  the  distance  from  the  center  of  the  earth  r  of  the 
current  position  P  of  the  rocket  is 

r-V(aV +  /?..)•+ (K  +  /?.,)*  +  (Z  +  /?..)'  (5.7) 


Obviously,  the  invariant  inequalities 
2  2  .01 
r  «R  are  all  established*  Because  of  this,  it  is  accurate  to 

make  use  of  the  height  h  =  R-R^^  by  which  the  rocket  leaves’  the 

surface  of  the  earth  in  the  inertial  coordinate  system 

*^o*o^o^o  handle  the  atmospheric  constants  / P ^ 

P/Po- 

The  expressions  (5.5)  and  (5.6)  for  the  velocity  V  and  the 

position  in  the  dynamic  coordinate  system  OXYZ,  arrived  at  from 

going  through  transformations  of  the  parameters  and  ^  in  the 

inertial  space  O^X^Y^Z^,  are  important  results.  They  clearly 

show  that,  in  O^X^Y^Z^,  solving  kinetic  equations  is  capable 

of  satisfying  various  types  of  requirements.  For  example,  during 

rocket  test  flights,  ground  optical  transit  theodolites,  radar,  or 

other  similar  tracking  and  measuring  equipment  requires  V  and  TJ)  to  act 

as  a  reference  basis.  In  these  two  types  of  coordinate  system,  the 

results  from  simultaneously  solving  rocket  kinetic  equations  clearly 

demonstrate  the  accuracy  of  (5.5)  and  (5.6),  verifying,  in  the 

coordinate  system  OXYZ,  that  solving  the  rocket  kinetic 

o  o  o  o 

equations  is  correct. 
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VI.  GEOGRAPHICAL  COORDINATE  COEFFICIENTS 


37 


As  far  as  the  geographical  latitude  and  the  geographical 

longitude  A  of  the  position  P  for  the  current  location  of  the  rocket 

at  the  instant  t  is  concerned,  it  is  also  possible  to  make  use  of  the 

position  R  in  O^X  to  accurately  calculate  it  out. 

o  o  o  o  . 

Refering  to  Fig.l,  the  geographical  latitude  <p  ^  can  be  obtained 
from 


singj.=£//? 

(6.1) 

In  this,  R  is  the  distance  of  the  point  P  from  the  center  of  the  earth 
E.  ^  is  the  vertical  height  of  the  point  P  off  the  plane  of  the 
equator.  The  geographical  latitude  is 

In  reality,-  this  only  has  meaning  on  the  s-urface  of  the  earth. 

However,  cp  ^  has  significance  everywhere.  The  line  connecting  the 
center  of  the  earth  E  and  P  intersects  with  the  surface  of  the  earth 
at  point  P' ,  In  such  a  case,  the  distance  of  the  point  P'  from  the 
center  of  the  earth  is 


y?'=  i?,  >/l  —  «• /V  1— c*'COsV» 


(6.3) 


As  far  as  the  rocket's  reaching  point  P  in  space  at  the  current 

instant  t  is  concerned,  in  terms  of  the  take  off  point  O  's  inertial 

o 

space  span  or  tensile  angle,  ( unclear )  voyage  or  course 

angle.  Using  the  numerical  product  of  the  two  vectors  from  the  points 
0^  and  P  pointing  toward  the  center  of  the  earth  E,  it  is  expressed 
as 


»Rty  Z%Rt, 

R,R 


fS 


(6.4) 


After  going  through  the  time  t  at  the  earth’s  prime  meridian,  in 

inertial  space,  it  turns  from  EX  '  to  EX  ,  that  is,  about  the  axis 

of  autorotation  of  the  earth,  EZ  turns  through  the  angle  to  t. 

©  © 

In  the  same  way,  OXYZ  also,  from  the  originally  duplicate  location  at 

O  X  Y  Z  ,  rotates  through  the  same  angle  CC^t.  From  Fig.l, 
oooo  ^  ^  e 

one  can  see  that  the  difference  in  longitude  of  the  point  P  relative 
to  point  0^  in  inertial  space  is 


(6.5) 


On  the  spherical  surface  triangle  EO^QP',  respectively  making  use  of 
cosine  theorems  and  sine  theorems,  it  is  possible  to  solve  for  an 
expression  for  geographical  longitude  A.  .  From  the  cosine  theorem  of 
the  course  or  angle  of  travel  unclear)' 

•  sin  cos  A/ 


Because  of  this,  one  obtains  the  formula  to  calculate  the  geographical 
longitude  A  rocket's  current  position  P 

cos (1 -!,+(»,/)  =  (cos/,  -  sin<p.,-  sin(p,)/(co6<p.,-cos<p.)  (6.7) 


In  this,  X  is  the  latitude  of  the  O  point  at  the  instant  of  the 
o  o 

rocket's  take  off,  t  =  0  sec.  In  order  to  accurately  determine  the 
quadrant  of  A  ,  on  the  same  spherical  surface  triangle  EO^QP',  one 
makes  use  of  the  cosine  theorem.  It  is  then  possible  to  solve  the 
expression  sin  (A“  +  CO^t).  From  Fig.l,  it  is  possible  to 

see  that  ^  ^  is  also  the  actual  geographical  longitude  of  tne  origin 
point  0  of  the  dynamic  coordinate  system  OXYZ.  The  geographical 
longitude  of  initiation  points  or  origins  and  0,  xc  ~  ^o 

In  this. 


xc 


/i«  *  sin  2B» 


(6.8) 


This  is  perpendicular  deviation.  However, 


is  the  rate  of  flattening  of  the  earth.  R  and  R.  are, 

®  J 

respectively,  the  equatorial  radius  and  the  polar  radius  of  the  earth. 
The  azimuth  angle  of  the  current  position  P  is 


y,=j4,+/i 


(6.10) 


In  this,  f is  the  angle  of  offset  or  deviation  relative  to  the  main 
plane  of  the  flight  foe  the  current  time  position  of  the 

rocket  P,  that  is,  the  included  angle  between  O^P  and 
Obviously, 


sin*'^- 


r) 


(6.11) 


From  this,  it  is  possible  to  obtain  a  precise  determinatinon  of  the 
quadrant  of  X  . 

CONCLUSIONS 


The  practical  realization  of  the  study  and  design  of  the  inertial 
guidance  systems  of  our  carrier  rockets  clearly  demonstrates  that,  in 
the  inertial  coordinate  system  which  is  illucidated 

by  this  article,  the  solutions  of  rocket  kinetics  and  the  derivations 
of  various  individual  formulae  are  correct.  This  is  particularly  true 
for  the  motion  parameters  V  and  y]  as  well  as  geographical  coordinates 
B  and  X  and  is  completely  in  line  with  the  results  obtained  in  the 
dynamic  coordinate  system  OXJfZ  from  solving  kinetic  equations. 

20 


38 


-IT. 

Making  use  of  the  inverse  matrix  B  =B  ,  taking  carrier 

rockets'  actual  flight  course  track  parameters  and  transforming  them 

into  O^X  y  Z  /  the  carrying  out  of  checks  and  tests  on  the 
o  o  o  o 

errors  of  inertial  instrumentation  and  the  anaFysis  of  other  errors  in 

guidance  systems  must,  of  necessity,  be  feasible. 

12  3  4 

There  are  a  number  of  representative  works  '  '  '  and  treatises 
which,  in  the  dynamic  coordinate  system  OXYZ,  which  rotates  following 
the  earth,  set  up  and  solve  kinetic  equations  of  rockets,  handle 
various  types  of  problems  associated  with  flight  force  mechanics, 
calculate  aerodynamic  forces  with  relative  ease,  make  observations  and 
measurements  as  well  as  tracking  rocket  motions  without  the  need  for 
transformations.  With  regard  to  the  extreme  importance  of  motions  and 
errors  associated  with  the  study  of  carrier  rocket  inertial  guidance 
systems  in  inertial  space  as  well  as  inertial  instrumentation,  the 
methods  studied  in  this  article  have  value.  In  particular ,  writing 
them  down  into  an  article  is  convenient  for  use  as  a  reference.  In 
conjunction  with  that,  they  have  been  even  more  deeply  discussed  with 
coworkers. 

In  the  process  of  the  studies  in  this  article,  we  were  aided  by 
Li  Zhentao, -Wang  Guangmin,' and  other ‘ similar  comrades.  We  have  opted 
for  the  use  of  several  formulae  associated  with  rocket  kinetic 
equations  set-  up  by  them  in  the  OXYZ  system.  We  wish  to  express  our 
sincere  thanks  to  theml 
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